Bright solitons in Bose-Fermi mixtures 
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We consider the formation of bright solitons in a mixture of Bose and Fermi degenerate gases 
confined in a three-dimensional elongated harmonic trap. The Bose and Fermi atoms are assumed 
to effectively attract each other whereas bosonic atoms repel each other. Strong enough attraction 
between bosonic and fermionic components can change the character of the interaction within the 
bosonic cloud from repulsive to attractive making thus possible the generation of bright solitons in 
the mixture. On the other hand, such structures might be in danger due to the collapse phenomenon 
existing in attractive gases. We show, however, that under some conditions (defined by the strength 
of the Bose-Fermi components attraction) the structures which neither spread nor collapse can be 
generated. For elongated enough traps the formation of solitons is possible even at the "natural" 
value of the mutual Bose-Fermi (*^Rb -""^K in our case) scattering length. 
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Quantum degenerate atomic gases offer unique oppor- 
tunities for investigating nonlinear properties of mat- 
ter waves. Atomic Bose-Einstein condensation, first ob- 
served in 1995 1], makes possible sources of sufficiently 
intense atomic matter waves allowing for the observation 
of nonlinear effects similar to those present in nonlin- 
ear optics. An important difference is that while the 
light must propagate through the nonlinear medium to 
produce the nonlinear effects, the nonlinearity is auto- 
matically present for atomic matter waves due to in- 
teractions between atoms. First confirmation of nonlin- 
ear atomic matter waves properties came in 1999 when 
the four-wave-mixing phenomenon was observed in the 
Bose-Einstein condensate of sodium atoms 0. The sec- 
ond was a series of experiments discovering the solitonic 
behavior of atomic condensates. First, it was demon- 
strated with the help of the phase imprinting technique 
that the condensate wave function for repulsive in- 
teractions exhibits the dark solitonic solutions Q. Next, 
also bright matter-wave solitons were created experimen- 
tally but this was achieved by other methods. In one of 
them, the interactions between atoms are changed from 
repulsive to attractive using the magnetically tuned Fes- 
hbach resonances technique [3, the other one makes 
use of possibility of changing the sign of effective mass 
in the periodic potential (i.e., in the optical lattice), a 
concept developed in Ref. 0. 

In Ref. we proposed a new way for generating 
bright solitons in repulsive Bose-Einstein condensates. 
The idea was to immerse the condensate in a degener- 
ate gas of fermionic atoms which attract bosons. It turns 
out that for strong enough attraction between bosons 
and fermions the system starts to behave as a mixture of 
two effectively attractive gases even though none of the 
components alone is an attractive system. Therefore the 
formation of bright solitons becomes possible. These soli- 
tons are two-component peaks with much smaller number 
of fermions than bosons. However, in a three-dimensional 
space the atomic gases with attractive forces might show 
an instability due to the collapse. This is the main goal 
of this paper to answer the question whether solitons in 



a Bose-Fermi mixture can be generated despite of the 
existing instability driven by the attractive forces. 

The system, we consider, is a Bose-Fermi mixture con- 
fined in a three-dimensional trap at zero temperature. 
We focus on the mixture of *^Rb and ''"K atoms in their 
doubly spin polarized states. At zero temperature this 
system is described by the many-body wave function 
^'(xi, ...jXAT^iyi, ...,yNp\t), where Nb {Np) is the num- 
ber of bosons (fermions). To find the time evolution of 
the system we impose an assumption according to which 
the many-body wave function is given in terms of single- 
particle orbitals in the following way 

Nb 

^(xi, ...,XAr^;yi, ...,yjv_p;i) = ]J <^^^^ (xj, i) 
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In fact, it is the simplest form that preserves the appro- 
priate symmetry conditions for both the bosonic part and 
the spin-polarized fermionic subsystem. 

To find the equations that govern the evolution of the 
single-particle wave functions we follow the formalism 
based on the Lagrangian density. Since the fermionic 
sample is spin-polarized the only interactions included 
are the contact boson-boson (assumed to be effectively 
repulsive) and boson-fermion (assumed to be attrac- 
tive) ones. The time-dependent many-body Schrodinger 
equation is fully equivalent to the following Lagrangian 
density 
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treated as a function of the many-body wave func- 
tion \E'(xi, ...,XjvB;yi, ■■■,yNFT't) its time and spatial 
derivatives. 

Therefore, the many-body wave function of the form 
is inserted into the Lagrangian density Next, 
we apply the mean-field approximation by integrating 
each term of over all but one particular bosonic or 
fermionic coordinate (when the interaction term is con- 
sidered both the one bosonic and one fermionic coordi- 
nates have to be left out). Utilizing the orthogonality 
properties of the fermionic orbitals (as well as neglecting 
1 in comparison with the number of bosons Nb) one gets 
the effective Lagrangian density 
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which depends on single-particle bosonic and fermionic 
orbitals and their spatial and time derivatives. Here, we 
already introduced the interaction strengths gs and gsF 
that determine relevant interatomic interactions and are 
related to the scattering lengths as and ubf via formulas 
A-Kfi^aB/mB and 2'Kh^aBF / IJ-, respectively (mg is a mass 
of *^Rb atom and /i is a reduced mass of ^^Rb and 
atoms). 

Having the effective single-particle Lagrangian density 
one can easily derive the corresponding Euler-Lagrange 
equations (here, j = 1, 2, Nf) 
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Neglecting the interaction between bosons and fermions 
results in a separation of the above equations; the first 



of them becomes then the Gross-Pitaevskii equation for 
a degenerate Bose gas whereas the rest describes the 
system of noninteracting fermions. However, when the 
Bose and Fermi components attract each other strongly 
enough the last term in the first equation can dominate 
over the term describing the bosonic mean-field energy 
causing effective attraction between bosons. We are go- 
ing to investigate further this possibility. 

We solve numerically the Eqs. Q for a mixture of 
^^Rb-^°K atoms confined in a three-dimensional elon- 
gated harmonic trap as well as in an atomic waveguide 
(obtained by releasing the axial confinement of the har- 
monic trap). We start from getting the ground state 
of a three-dimensional mixture in a harmonic trap. To 
this end, we pick up the atomic orbitals at zero time as 
they were the stationary solutions of Eqs. I@J assuming 
no interaction between bosons and fermions {asp = 0). 
Therefore, the bosonic wave function is the solution of 
the Gross-Pitaevskii equation, whereas the fermionic or- 
bitals are successive (with respect to the energy) eigen- 
states of the harmonic potential. Because of numerical 
limitations, the number of fermionic atoms is small (at 
most 15 in our case) and since the trap is elongated only 
states with higher axial quantum number are populated 
(the radial quantum number equals zero). Next, we adia- 
batically increase the attraction between Bose and Fermi 
components. 
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FIG. 1: Axial density profiles of a three-dimensional Bose- 
Fermi mixture consisting of 100 bosons (^'^Rb atoms) and 10 
fermions (*°K atoms) confined in the harmonic trap with axial 
frequency 2-kxI Hz and the radial one 2% x 300 Hz for bosonic 
component (corresponding trap frequencies for fermions are 
multiplied by a factor 1.473). Interactions between bosons are 
determined by the "natural" scattering length as = 5.8 nm 
whereas the strength of *^Rb atoms -*''K atoms attraction 
is being changed. Dark (bright) shaded areas correspond to 
fermions (bosons). At zero attraction between fermions and 
bosons, fermions form a broad density distribution due to 
the Pauli exclusion principle. For strong enough attraction, 
however, fermions are pulled inside the bosonic cloud. 

In Fig. ^we plot the axial density (in units if l/a^^, 
where aho = (S./tof'^z^'')"'^^^) of the ®^Rb-^"K mixture 
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with 100 bosons and 10 fermions confined in the har- 
monic trap. The densities are normahzed to one. The 
trap frequencies for bosons are 27r x 1 Hz {2tt x 300 
Hz) in axial (radial) direction (corresponding frequen- 
cies for fermions are multiplied by a factor 1.473). For 
weaker attraction between bosonic and fermionic atoms 
fermionic cloud broadens due to the Pauli exclusion prin- 
ciple. Stronger attraction, however, can overcome the 
repulsion caused by the Pauli exclusion principle and all 
fermions are pulled inside the bosonic cloud. Increasing 
further the attraction results in an extremely fast growth 
of the density at the center of the trap. In other words, 
the collapse of the mixture happens. It is shown in Fig. 
13 This phenomenon was already observed experimen- 
tally 9] as well as extensively studied theoretically T^. 
It is well known that such a behavior is attributed to 
the three-dimensional realm and does not happen in a 
one-dimensional geometry. 



150 

DENSITY 

y5 (osc. units) 




-10 -20 -30 -40 -50 -60 
Sbf (nm) 



FIG. 2: Illustration of a collapse phenomenon. All parameters 
but the number of fermions (that equals 15 now) are the same 
as in Fig. Q The lower panel shows the bosonic (solid line) 
and fermionic (dashed line) densities at the center of the trap 
as a function of the scattering length asp- The onset of the 
collapse at about asp = — 60nm is clearly visible. 

Having the ground state of the ^^Rb-'*°K mixture in 
the harmonic trap we can now reload the system to the 
atomic waveguide just by opening the trap in axial di- 
rection. Results of this procedure are presented in Figs. 
13 and |31 The response of the system to the releasing of 
the axial confinement depends on the value of the scat- 
tering length asF- It turns out that for a particular trap 



geometry there is an interval of values of qb f with such 
a property that ii ubf lies within this interval the two- 
component single-peak structure is formed which does 
not spread nor collapse. This case is shown in Fig. |21 
where the axial confinement is linearly open in approxi- 
mately 0.5 s. The mixture looses 6 bosons while both den- 
sities fit to each other and then apparently stays in the 
center of the trap without changing its shape. To check 
whether the mixture forms a soliton (or at least a solitary 
wave) we forced it to move by pushing along the waveg- 
uide direction the bosonic component. It is clear from 
Fig. 01 that bosons are able to pull fermions and both 
components stay together during the movement. Thus 
the mixture behaves as a solitary wave although some 
breathing of the structure is observed. The axial shapes 
of the densities resemble rather step-like function than 
the familiar solitonic secans-hyperbolicus shape. We find 
this shape as a typical (i.e., present in all geometries in- 
vestigated) while changing the scattering length a^f by 
approaching the edge of the window supporting the exis- 
tence of solitons, corresponding to the weaker attraction. 
For stronger attraction the shape is "more" solitonic (see 
Fig. EJl. The soliton in Fig. |21 is 40 /im long and the 
bosonic density is of the order of 10^^ cm"'^. 
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FIG. 3: Two-component soliton standing at the center of the 
trap. The mutual scattering length equals asF = — 40nm. 
The axial confinement (the radial one is still kept on) is re- 
moved linearly within 0.5 s. Bosonic and fermionic densities 
fit to each other in approximately Is (the system looses about 
6 bosons within this time) and later on the single peak two- 
component soliton is formed. Radial trap frequency is the 
same as that in the case of Fig. Q 

When the value of the scattering length asF lies out- 
side the window just discussed (i.e., the attraction be- 
tween bosons and fermions is weak enough) we see the 
spreading of bosonic and fermionic clouds (Fig. [SJ. Here, 
ttBF = — 30nm and the attraction is too weak to glue 
the bosonic and fermionic clouds. Both clouds spread 
and eventually break into several smaller droplets. Fig. 
|S1 shows that at least outermost droplets (consisting of 5 
fermions and 30 bosons) demonstrate the solitonic behav- 
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FIG. 4: Illustration of the solitonic character of a structure 
obtained by releasing the axial confinement as in the case of 
Fig. 121 The structure is pushed along the axial direction 
(direction of no confinement) and is moving maintaining its 
shape (except small "breathing"). This is the case when the 
scattering length ubf lies within the window supporting the 
existence of bright solitons. 

ior moving without changing their shapes. This observa- 
tion supports the view that weaker attraction allows only 
for two-component solitons with small number of atoms 
(as opposed to what is reported in Ref. flU). 
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FIG. 5: Nonsolitonic behavior of bosonic and fermionic 
clouds. The attraction between bosons and fermions is not 
strong enough {olbf ~ —30 nm) and both densities are spread- 
ing. Later on, however, the clouds break into several pieces 
that move further without changing their shapes illustrating 
the rule that at weaker attraction solitons with smaller num- 
ber of atoms are allowed. Here, the initial number of fermions 
equals 15. 

Finally, we consider the collision of two bright two- 
component solitons within the waveguide (Fig. EJ. Both 
of them contains 6 fermions and 100 bosons. The right 
one stays at rest whereas the left one is forced to move by 
pushing the bosonic component. The left soliton is mov- 
ing with the velocity 0.12 mm/s. Moreover, the phase tt 



is added to the bosonic part of the soliton at rest, other- 
wise solitons merge when they meet at the center of the 
waveguide and collapse. This kind of behavior is simi- 
lar to what happens for pure bosonic solitons governed 
by the Gross-Pitaevskii equation. In that case solitons 
repel or attract each other depending on the initial rel- 
ative phase (l^ . If this relative phase equals tt the pure 
bosonic solitons always repel each other. This is what we 
observe in our numerical simulations, an example is given 
in Fig. El After collision the right soliton starts to move 
whereas the left one remains almost at rest. Usually, the 
transfer of some mass happens during the collision. 
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FIG. 6: Collision of two solitons, each of them consisting 
of 6 fermions and 100 bosons. The left soliton is moving 
towards the right one being initially at rest. After collision 
(accompanied by the transfer of some atoms) the right soliton 
is moving and the left one is staying almost at rest. In this 
case ubf = — 120nm and the trap parameters are 27r x 1 Hz 
in the axial direction and the 2tt x 30 Hz for the radial one 
for bosonic component. 

To summarize our numerical work we present in Table 
^ some numbers characterizing the formation of single- 
peak two-component bright solitons. In particular, two 
last columns give approximate values of the critical mu- 
tual scattering lengths a^p (the mixture spreads for the 
weaker attraction than a^p) and the values of a'^p de- 
termining the onset of the collapse for various radial con- 
finements. Thus, there exist windows supporting the ex- 
istence of stable structures that maintain their shapes 
while moving and may survive the collisions (although 
some atoms are usually transfered during the collision). 
Such structures can be called solitons. The critical value 
of the scattering length can be also estimated as in the 
Ref. 0. Based on the first equation of a set of Eqs. 
^ the condition for the system to enter the phase when 
the bosonic component becomes the gas of effectively at- 
tractive atoms (and hence allowing for the generation of 
bright solitons) is written as gsns = \9'bf\''^f, where 
?i_B = iVslyj^-^^p and np = J^i'" I'/'i^'P ^-^^ the bosonic 
and fermionic densities at the center. This criterion gives 
the values: — 86nm, — 40nm, and — 14nm for the radial 
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confinements: 27r x 30 Hz, 27r x 300 Hz, and 27r x 1000 Hz, 
respectively, what remains in a good agreement with the 
values obtained directly from the numerical simulations. 



TABLE I: Parameters characterizing the formation of single- 
peak two-component bright solitons. 



i/P'pz] 


Nb 


Nf 






30 


100 


6 


-90 


-135 


300 


94 


10 


-35 


-62 


1000 


64 


15 


-17 


-44 



In Table we put also the numbers of fcrmionic and 
bosonic atoms within the single-peak solitons. They 
are obtained in the following way. For each geometry 
(i.e., the radial confinement) we start from the three- 
dimensional harmonic trap ground state consisting of 100 
bosons and 10 or 15 fermions. After releasing the ax- 
ial confinement the excess of atoms (bosons or fermions) 
flows out of the system. For example, for the radial fre- 
quency Vr — 30 Hz the mixture looses 4 fermions and 
then behaves as a soliton (to get more fermions within 
the peaks the initial number of bosons would have to 
be bigger). For stronger radial confinement the mix- 
ture expels the excessive bosons. We have checked that 
for ' = 2000 Hz there exists a soliton at the natural 
value of the scattering length aBP (equal to — 13.2nm, 
[T3|). Our simulations show that at small values of \aBF\ 
only solitons with small numbers of atoms exist, the ex- 



cessive bosons or fermions just leak out of the system. 
Our preliminary results based on the three-dimensional 
variational analysis within the hydrodynamic model of a 
Bose-Fermi mixture Tij also demonstrate that solitons 
with large numbers of atoms require stronger attraction 
between fcrmionic and bosonic components. 

In conclusion, we have shown that bright solitons can 
be generated in a Bose-Fermi mixture trapped in a three- 
dimensional elongated harmonic potential. For that the 
attraction between bosonic and fcrmionic atoms has to 
be strong enough changing in this way the repulsive in- 
teractions among bosons to the attractive ones. Simul- 
taneously, the attractive forces induced in the bosonic 
cloud has to be weak enough to avoid the collapse. It 
turns out that there exist windows in the values of the 
rubidium-potassium scattering length within which the 
single-peak two-component structures neither spreading 
nor collapsing (therefore called by us solitons) can be 
formed. Those windows are shifting towards the region 
of less negative Bose-Fermi scattering length (weaker at- 
traction) while elongating the trap making thus possible 
the formation of solitons even at the natural strength of 
the boson-fermion attraction. 
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